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Abstract 

We describe the LDA bandstructure of YBa2Cu3 07 in the e_p ± 2 eV range using orbital 
projections and compare with YBa2Cu4 08. Then, the high-energy and chain-related degrees 
of freedom are integrated out and we arrive at two, nearest-neighbor, orthogonal, two-center, 
8-band Hamiltonians, and , for respectively the even and odd bands of the bi-layer. 
Of the 8 orbitals, Cu^2_j,2, 02^, 03j,, and Cus have a character and Cu^z, Cu^^, 02z, and 
03^ have tt character. The roles of the Cus orbital, which has some Cu3^2_i character, and 
the four tt orbitals are as follows: Cus provides 2nd- and 3rd-nearest-neighbor (t' and t") 
intra-plane hopping, as well as hopping between planes {t±). The 7r-orbitals are responsible 
for bifurcation of the saddle-points for dimpled planes. The 4-(7-band Hamiltonian is generic 
for flat Cu02 planes and we use it for analytical studies. The k||-dependence is expressed 
as one on m = (cosbky + cos akx) /2 and one on v = (cosbky — cos akx) /2. The latter arises 
solely through the influence of Cu^. The reduction of the cr-Hamiltonian to 3- and f-band 
Hamiltonians is explicitly discussed and we point out that, in addition to the hoppings com- 
monly included in many-body calculations, the 3-band Hamiltonian should include hopping 
between all 2nd-nearest-neighbor oxygens and that the f-band Hamiltonian should include 
3rd-nearest-neighbor hoppings. We calculate the single-particle hopping between the planes 
of a bi-layer and show that it is generically: t± (ky) f» 0.25 eY-v'^ (1 - 2ut' /t) ^ . The hop- 
ping through insulating spacers such as (BaO)Hg(BaO) is an order of magnitude smaller, 
but seems to have the same k||-dependence. We show that the inclusion of t' is crucial for 
understanding ARPES for the anti- ferromagnetic insulator Sr2Cu02Cl2. Finally, we estimate 
the value of the inter-plane exchange constant Jj_ for an un-doped bi-layer in mean-field 
theory using different single-particle Hamiltonians, the LDA for YBa2Cu306, the eight- and 
four-band Hamiltonians, as well as an analytical calculation for the latter. We conclude that 
Jj. 20 meV. 

1 Introduction. 

For the HTSC materials practitioners of the LDA now basically agree about their results. The 
major uncertainties come from the compositional and structural data rather than from the compu- 
tational techniques. For the various materials, the energy bands originating from the Cu02-planes 
are fairly similar, merely the position of the Fermi level changes. 

For stoichiometric HTSC's such as YBa2Cu307 the LDA has been successful in reproducing or 
predicting normal-state structural and phononic properties [1, 2, 3], optical excitation spectra at 
high energies (Tiw > 1 eV for metals, Tiw > 3 eV for Mott-Hubbard insulators) [4], and fc^-averaged 
Fermi surfaces e (k||)=0 for the metals [5, 3]. The most surprising LDA prediction for HTSC's 
with dimpled Cu02 planes (YBa2Cu307 [1, 3], YBa2Cu306 5 [6]) has been that of bifurcated 
saddle-points close to the Fermi level. Extended saddle-points near ep were subsequently found in 
ARPES [7], but whether these are related to the bifurcated LDA saddle-points through "Fermion 
condensation" [8], or are caused by anti-ferromagnetic spin-fluctuations [9, 10], or are due to 
surface effects [11], remains to be seen. The LDA has finally provided reasonable estimates of the 
electron-electron (e.g. Ucu d ~ 7 eV) [12, 13, 14] and electron-phonon interactions (e.g. Xg-ph ~ 1 
and softening of the Raman-active dimpling modes for T < Tc) [1, 15, 3, 16]. Xg-ph ~ 1 might 
account for < 40 K, but hardly for 100 K. 

However, the LDA band structures for HTSC's are complicated, rarely interpreted, not accu- 
rate - or even relevant - below 50 meV, and are not delivered in a form useful as the single-particle 



term of a correlated model Hamiltonian which describes the low-energy excitations. As a result, 
most theorists neglect the LDA band structure, or at least its non-trivial details, and use the 
simplest possible Cu-Cu one-band, two-center, orthogonal, tight-binding model with hopping in- 
tegrals t and t' between only nearest and next-nearest Cu neighbors (see Eq. (16)). A few brave 
theorists [10, 9] use the Emery three-band model. 

We have therefore integrated out of the LDA band structure for YBa2Cu307 the high-energy 
and chain-related degrees of freedom. In a first step we used the LMTO downfolding technique [17] 
to derive a two-dimensional, nearest-neighbor, two-center, orthogonal, tight-binding Hamiltonian 
with eight orbitals per Cu02 plane [18]. Of these eight Cu02 orbitals, four (Cuj;2_j,2, 02^, 03y, 
and Cus) have a character and four (Cu^;^, Cuj,^, 02^, and 03^) have tt character with respect to 
the plane. Comparisons with full LDA bands calculated for other structurally well-characterized 
HTSC materials convinced us, that this eight-band Hamiltonian is generic for LDA plane bands. 

The present paper starts with a qualitative description of the full LDA band structure for 
YBa2Cu307 and YBa2Cu408 in the ep ±2 eV region and emphasizes the c/jazw-related features. 
Then we proceed with a quantitative description of the plane bands by means of the eight-band 
Hamiltonian. The roles of the Cuj orbital, which possesses some Cu3^2_]^ character, and the 
four TT-orbitals are specifically explained: Cuj provides 2nd- and 3rd-nearest-neighbor (t' and t") 
intra-plane hopping, as well as hopping between planes (t_i). The 7r-orbitals are responsible for 
bifurcation of the saddle-points. We explicitly show how further high-energy degrees of freedom 
may be integrated out. The four cr-bands are studied analytically. Their ky -dependence is most 
simply expressed as one on m = (cos hky + cos ak^) /2 and one on t; = (cos hky — cos ak^) /2, where 
the latter arises solely through the infiuence of Cuj . The reduction of the cr-band Hamiltonian to 
the commonly used 3- and 1-band Hamiltonians is explicitly discussed, and we point out that, to 
be consistent with the LDA bands, the former should include hopping between all 2nd-nearest- 
neighbor oxygens and the latter should include 3rd-nearest-neighbor hopping. As a first simple 
application, we calculate the single-particle hopping between the planes of a bilayer and show that 
it is generically t±_ (ky) ~ 0.25 eV-t;^ (1 — 2ut' /t) ^ . The hopping through insulating spacers such 
as (BaO)Hg(BaO) is an order of magnitude smaller, but seems to have the same ky-dependence. 
These hopping integrals are relevant for c-axis transport and their squares enter the Inter-layer 
Pair- Tunnelling model [19]. Then we show that the inclusion of t' is crucial for understanding 
ARPES [20] for the anti-ferromagnetic insulator Sr2Cu02Cl2. Finally, we estimate the value of 
the inter-plane exchange constant Jj_ for a non-doped bi-layer in mean-field theory using different 
single-particle Hamiltonians, the LDA for YBa2Cu306, the eight- and four-band Hamiltonians, as 
well as an analytical evaluation for the latter. We conclude that Jj_ ~ —20 meV. 

2 LDA band structures of YBa2Cu307 and YBa2Cu408. 

The full LDA band structure for YBa2Cu307 [3] between the ky points (0,0), (f,0) , , f ) , 
(0,f) , (0,0), and (f ,f) is shown in Fig. 1 for k,={) (FXSYFS), and in Fig. 2 for k, = ^ 
(ZURTZR). Also shown are the cross-sections of the Fermi surface with the A;^=0 and planes. 
Each of the small subfigures gives the projection onto a particular O p or Cu d orbital: The radius 
of the circle around the e„ (ky) -point is proportional to the weight of that orbital in ipn (ky). 
Note that Cu and O orbitals have different normalizations in the figures so that weights cannot 
be compared between copper and oxygen. 03 are the plane-oxygens running along h, which is the 
direction of the chain. 02 are plane-oxygens along a, Cu2 are plane coppers, 04 apical- and 01 
chain-oxygens, and Cul are chain-coppers. 

The FS has four sheets denoted a, b, c, and s. The essential wave functions for the a and 
b sheets are respectively the antibonding (odd, -) and bonding (even, -|-) linear combinations 
for the bilayer of the plane pda anti-bonding orbitals (02j;-Cu2j;2_j,2-03j,). Essential for the c 
sheet is the chain pda anti-bonding orbital made from Olj, and the pda anti-bonding dumbbell 
orbital (04^-Cul^2_j,2-04^). Note that the orbitals denoted in the figures are d^z^-i and that 
dz2_y2 = {\/?>d2,z2-i + dj;2_y2^ /2. Tlic orbitals essential for the s sheet, the small stick pocket 
along SR, are two out of the three chain pdir orbitals, which are the two pdir anti-bonding dumbbell 



orbitals (04j;-Culj;^-04j; and 04j,-Culj,^-04j,) plus the anti-bonding linear combination of 01^ 
with the pdiT bonding dumbbell orbital 04j,-Culj,^-04j,. The anti-bonding linear combination of 
the two latter orbitals, which both involve OAy, is the essential wave function for the stick. 

If the one-electron wave functions stay coherent from chain to chain, as it is assumed in the 
LDA for a perfect crystal, then for A;^=0 there is neither hybridization between the a and c sheets, 
nor between the b and s sheets. Conversely, for ^^ = 7 there is neither hybridization between the 
b and c sheets nor between the a and s sheets. For a general k;,, all chain-plane hybridizations 
are allowed and give rise to a fc^-dispersion of several hundred meV's. 

The Fermi level is approximately 20 meV above the odd plane-band saddle-points, which are 
bifurcated away from X and Y to respectively Xi^FX and Y±^FY. Since b/a 1.015 and 
since the ratio between the Cu2j;2_j,2-03j, and Cu2j;2_j,2-02j; hopping integrals scales roughly 
as {h/ a) ^ , we expect the saddle-point near X to lie more than 100 meV higher than the one 
near Y. However, the calculated energy difference is only 10 meV. The reason is that the odd 
plane band is repelled along FY by the below-lying chain pdi: stick band {OAy and 01^). Further- 
more, the hybridization with the chain pda band for k^, ^ Q makes the plane-band saddle-points 
three-dimensional: The saddle-point related to (X ± ;|FX, -15 meV) is (U, -|-80 meV) and the one 
related to (Y ± ^FY, -25 meV) is (T ± ^FT, -300 meV) . From the fi gures we see that not only 
the odd, but also the even plane band has bifurcated saddle-points, but these are about 0.5 eV 
below ep- 

Whether the above-mentioned fc^-dispersion persists in real YBa2Cu307 is unclear: In the 
douhle-cYiwi compound YBa2Cu408, the FDA dispersion is an order of magnitude smaller than 
in YBa2Cu307 because single-electron hopping through a double chain is frustrated by the fact 
that 01 in one chain plays the role of 04 for the other, and yet, the physical properties of the two 
compounds are quite similar. 

In YBa2Cu408 the FDA saddle-points of the odd plane band are (Xi^FX, -20 meV), (U, 
meV), (Yzb^FY, -160 meV) and (T, -180 meV). The difference between the energies near X and 
Y is thus as expected when the chain-plane hybridization is small and h/a 1.013. The even-odd 
splitting of the plane bands caused by hopping across the double-layer is ~0.5 eV near X or Y, 
exactly as in YBa2Cu307. 

3 8-band Hamiltonian. 

The chain degrees of freedom have low energy, but they are not generic because not all HTCS's 
have chains. In order to derive a generic low-energy few-band Hamiltonian for the planes, we have 
therefore downfolded the FDA multi-band Hamiltonian for YBa2Cu307 and A;^=0 to a two-center, 
orthogonal, tight-binding Hamiltonian iJ^ (ky) with eight odd bilayer orbitals interacting through 
nearest-neighbor hopping. Similarly, we have downfolded the FDA Hamiltonian for fc^ = f- to an 
eight-band Hamiltonian iJg" (ky) for the even bilayer orbitals [18]. Over a range of ±1 eV, 
and Hg accurately describe the odd and even plane bands, respectively, as they are seen in Figs. 
1 and 2. The weak hybridization with the chain pdi: bands has not been removed from, but is 
folded into these 8-band Hamiltonians. 

The downfolding consisted of first integrating out all degrees of freedom with very high energies 
(e.g. Ba spdf , Y spdf , Cu p, O sd). This was done using the FMTO downfolding procedure (see 
also Sect. 4). Thereafter, all remaining orbitals other than the eight odd (even) plane oxygen and 
copper orbitals specified in the Introduction were deleted from the downfolded FMTO basis, and 
the odd (even) plane bands were recalculated. These eight-FMTO bands accurately reproduced 
the corresponding full FDA plane bands over a ±2 eV range. Finally, the non-orthogonality and 
long range of the eight odd (even) orbitals were transformed away by fitting the bands in the ±1 
eV range to an orthogonal, two-center, nearest-neighbor, tight-binding Hamiltonian which then, 
is the one desired. 

The reason for keeping in the Hamiltonian particularly those eight orbitals, is that these orbitals 
are the ones which after orthogonalization describe the FDA plane-bands accurately over a ±1 



eV range with nearest-neighbor hoppings only, that is, with the minimal number of parameters. 
The 8-band Hamiltonian is thus the one which is "chemically" meaningful and sufficiently simple 
that further degrees of freedom may be integrated out analytically, as we shall demonstrate in the 
following section. 

By transforming from orbitals even and odd with respect to the yttrium mid-plane of the 
bilayer to orbitals centered on a specific plane, one obtains a single-plane Hamiltonian, Hg (ky) , 
plus integrals for hopping between the two planes of the bilayer. 

Fig. 3 now synthesizes the eight plane-bands and gives the most significant parameters of Hg 
(tetragonal averages are for instance taken). Ref. [18] provides the accurate values of all the 
matrix elements, as well as the detailed expressions, and the comparison with the full LDA bands. 
The inter-plane hopping will be discussed separately in Sect. 8. 

Shown in the bottom left-hand panel of Fig. 3 are the four cr-orbitals looked upon from above 
the plane, \y) = 03y, \d) = Cu2j;2_j,2, \x) = 02^, and \s) = Cu2s. Below, we see, from the edge 
of the plane, the \d) orbital and two of the 7r-orbitals, \z) = 02^ and \xz) = Cu2xz- The orbital 
energies are given at the relevant points of the band structure shown in the left-hand panel, in eV 
and with respect to the energy of the Cu2j;2_j,2 orbital. Hence, £x=£y = £p is 0.9 eV below e^, and 
fs is 6.5 eV above. Moreover, eo2^=£03^ = is 0.4 eV above e^, and exz=^yz is 1 eV below. If we 
now, as shown in the figures at the bottom, include the pda hopping integrals txd=iyd = ipd = 1.6 
eV and the pdir hopping integrals tz^xz=iz,yz = 0.7 eV, we obtain the band structure shown in 
the panel above: The cr-orbitals give rise to a bonding, a non-bonding, and an anti-bonding 02^;- 
Cuj;2_j,2-03j, band plus a Cuj level (full lines), and the 7r-orbitals give rise to two decupled pairs 
of bonding anti-bonding bands which disperse in either the x oi y direction (stippled lines). The 
anti-bonding pda band, which will develop into the conduction band, has saddle-points at X (and 
Y) which are well above the top of the 7r-bands, and which are isotropic in the sense, that the 
absolute values of the band masses in the x and y directions are equal. This means that the FS 
at half-filling is a square with corners at X and Y. 

In the middle panel, the a and tt bands remain decupled, but we have introduced the G\ig-02x 
and Cus-03j, hoppings {tgp = 2.3 eV), as well as the tiny 03^-02^ hopping {tzz = 0.06 eV). The 
latter is the only one reaching beyond nearest neighbors and it merely lifts a degeneracy of the 
anti-bonding 7r-bands at F. The strong coupling of the remote Cuj orbital to the pda orbitals has 
the pronounced effect of depressing the conduction band near the saddle-points at X and Y, and 
thereby increasing the mass towards F and decreasing it towards S (tt, tt) . The FS passing through 
X and Y will therefore correspond to finite hole-doping and will bulge towards F (see the left-hand 
panel of Fig. 7). That the Cuj orbital, which has the azimuthal quantum number m^=0, may mix 
with the conduction band at X (and Y), but not at ^FS (f-, f-) , can be seen from the symmetries 
of the corresponding two anti-bonding pda wave functions shown schematically in Fig. 4. That 
there is, in fact, substantial by-mixing of Cuj character in the full LDA wave function at X, 
but none at ^FS, is proved by the \il> (r)p-contours shown in respectively the left and right-hand 
parts of Fig. 5. The by-mixing of Gu^z^-i character (which also has m^=0-symmetry) to the 
so-called Cuj orbital, is seen to be small. The net downwards shift of the saddle-point energy is 
the result of down-pushing by pure Cuj and weak up-pushing by pure Gu^z^-i, whose energy is 
below the saddle-point. In Sect. 6 we shall give analytical descriptions of the conduction band in 
this dpa-spa approximation. 

As a result of the hybridization with the Cuj orbital, the saddle-points of the anti-bonding 
(T-band straddle off the top of the appropriate 7r-band, so that even a weak dimple or buckle of 
the planes will introduce noticeable hybridization between the a and tt bands. This is seen in the 
right-hand panel where we have turned on the weak Cuj;2_j,2-02^ and Cuj;2_j,2-03^ hoppings 
(tzd = 0.24 eV cx sin (5). These couplings become allowed when there is a finite angle 6 7° 
in YBa2Cu307 and YBa2Cu408) between the Cu-0 bond and the plane of the 2D-translations. 
Since the anti-bonding pda and Cuj;^-02^ pdir orbitals can only hybridize between F and X, but 
not at X, the a-ir hybridization makes the saddle-point bifurcate away from X. 

In conclusion, the chemical eight-band Hamiltonian Hg (ky) , whose eigenvalues give the one- 



electron bands for a single plane, has the cr-block: 
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and the block mixing the a and tt orbitals for dimpled planes: 
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4 Integrating out high-energy degrees of freedom. 

In deriving the 8-band Hamiltonian from the LDA we have numerically integrated out many high- 
energy degrees of freedom. In the following, we shall integrate out further high-energy degrees of 
freedom starting from Hg or {a\ H \<t) . We now explain the Lowdin procedure: 

Partitioning an orthonormal basis into \i) and \j) , and down-folding the iJ-'-'-block yields: 

Hu'{e) = Hu' - Y^Hij [(Fii-e)-'l H^i, . (2) 

^ L J ii' 

J J 

This is exact in the sense that det [H — e] and det [H" (e) — e] have identical zeroes, namely 
the eigenvalues of H. If we now Taylor-expand the down-folded Hamiltonian about ep: H" (e) = 
H" (£f) + (£ — ep) H" (£f) + ..., we see that the energy dependence of the down-folded Hamiltonian 
to linear order arises from non-orthogonality of the modified i-basis, the overlap matrix being 
1 — H" (ef) • The higher-order energy dependence of H" (e) originates from energy dependence 
of the basis. The latter, and sometimes even the former energy dependencies can be neglected 
if the degrees of freedom to be integrated out have high energy, that is, if the energy range of 
interest lies far away from the eigenvalues of H^^ . An orthonormal low-energy Hamiltonian 7i 
may be obtained by orthogonalization of the modified i-basis, whose energy dependence has been 
neglected, e.g. 

n = eF+[l- H" (ep)] [H" (e^ ) - ] [l - H" (e^ )] . (3) 

5 3-band Hamiltonian. 

A popular tight-binding Hamiltonian for fiat planes is the 3-band Emery model. This is like 
((t| H |(t) in (1), but without the Cuj orbital, with a renormalized oxygen energy (cp e^) , and 
with 2nd-nearest-neighbor 02j;-03j, hopping. 



Projecting the Cuj orbital out of [<t\ H \<t) by means of (2) reveals that the values of and tpd 
are unchanged, but that 2nd-nearest-neighbor 02j;-03j, as well as 3rd-nearest-neighbor 02x-02x 
and 03j,-03j, hoppings, all of size tpp = t^^/ (fj — e) ~ 1.1 eV, must be added. Moreover, fp must 
be renormalized to = — 2tpp ~ — 3.0 eV which is, in fact, the standard value. All of this 
comes about, because the Cuj character (see Fig. 5) is now built into the tails of the neighboring 
02x and 03j, orbitals. 

According to the LDA, the Emery model should thus be modified to 
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Here we have neglected the non-orthogonality and, for the numerical values, we have taken the 
expansion energy, ep = e (X) + 1.53 eV. As long as e — e^? ^ — e^? ~ 5 eV, this is a good 
approximation. Orthogonalization by means of (3) would have introduced longer-range hoppings. 



6 The (7-band. 

In order to derive accurate and transparent expressions for the fiat-plane cr-bands e„ (ky) , or 
their constant-energy surfaces, we return to the 4-band Hamiltonian {a\H\a) in (1) and first 
Lowdin down-fold (2) the oxygen orbitals. The result is the following energy-dependent spa-dpa 
Hamiltonian: 
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where we have defined the new h^^-coordmates: 



(5) 



M = — (cos aky + cos akx) 



and V = — (cos aky — cos ak^) ■ 



(6) 



both limited to the range from -1 to and sketched in the left-hand side of Fig. 6. The origin 
of the Mt;-system is at ^FS and corresponding values of (k^, ky) and [u, v] are: ^FS (^, ^) [0, 0] , 
X , 0) [0, 1] , F (0, 0) [1, 0] , and S , f ) [-1, 0] . Readers used to think about pairing and gap 
symmetries, will recognize u as s- and v as rf-wave symmetry. The relation to the ky-coordinates 
X and y used in Ref. [18] is: I — u = x + y and v = x — y. 

The S(i-Hamiltonian (5) highlights the ky -dependence of the FDA cr-bands for a fiat Cu02 
plane as follows: Without s-mixing, the bands depend only on u (the saddle-points are isotropic, 
the FS for half filling is the square: m = 0, etc.) and the mixing element is proportional to v. 

Further Fowdin down-folding, this time of the spcr-block, yields the following energy-dependent 
1-band Hamiltonian: 
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(7) 



which is still equivalent with [<t\ H \<t) . When there is no mixing with Cu s, that is when tgp = 0, 
s{e) ^ 00 so that the term proportional to vanishes. The quadratic spcr-scattering function 



s (e) = (fs - e) (e - Cp)/ (2tspy 



(8) 



is fairly constant in the neighborhood of e = (e^ -|- £p) /2 2.8 eV, where it takes its maximum 
value [(fs — Cp) /4tsp]^ 0.65. For instance is s(1.5 eV) 0.57. 



The most elegant way of expressing the eigenvalues oi [<t \ H \<t) is as the roots of the secular 
equation, 
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because from its definition, 
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(9) 



(10) 

is a quadratic function of energy which depends only on the dpa scattering. For (9) to vanish, 
we see that d{e) must increase from to 2 in the region of the conduction band because s(e) 
is positive and because the bottom of the band is at F and the top is at S. For a given energy, 
e, we calculate the values of s(e) and d{e) from (8) and (10), whereafter the requirement that 
(9) be zero gives the exact constant-energy surface. For a given ky , on the other hand, we can 
not give an exact, explicit expression for the four eigenvalues, e„ but near the middle of the 

conduction band we may exploit the near-constancy of s (e) and are then led by (10) and (9) to 
solve the following equation: 
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For ease of notation, we have dropped the subscript || on kii 



(11) 



iteratively [i.e. e (k) e] 
shall do so from now on. 

With the parameter-values from Fig. 3, we easily obtain: e {X) = e (m=0, t;=l) 1.53 eV and 
e (^FS) = e(M=0,i'=0) 2.78 eV. The Cuj level thus pushes the saddle-point down by 1.25 eV. 
The energy of the bottom of the conduction band, e (F) = e{u=l,v=Q) = max(ep,e(j) = 0, and of 
its top, e (S) = e{u=-l,v=Q) 4.10 eV, are independent of s. 

The CU(j, Cus, and Op (02-1-03) characters of the conduction-band wave function are most 
easily obtained as the respective derivatives (9e (k) /Se^, (9e (k) /dep, and (9e (k) /deg. This follows 
from Ist-order perturbation theory. By differentiation of (9) and making use of 



rf(e) = 2i^%±^ and . (e) = 2 i^^ + ^il^^. 
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we obtain the simple results: 
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and \cp (e, k)p = 1 — |cs (e, k)|^ — |C(j (e, k)p . The Cu^ character is seen to be ~ 5 [1 ~ ••] ^ when the 
energy is far away from £p ~ e^. Moreover, the variation of Cu d character along a constant-energy 
surface, e.g. the FS, becomes negligible when 
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£-2.8 eV 



e-F 0.45 eV 



(15) 



and vanishes for the energy midway between the Op and Cuj energies. With our parameter-values 
and for e=1.53 eV, s/d 0.31, whereas for e=2.78, s/d = 0. In the following, we shall often 
neglect s, and when this can be done, the Cu^ character is independent of k and slightly larger 
than 0.5. The Cuj character is proportional to and, with our values of the hopping integrals 
and our s-value, \cg (e, k)p ~ 0.12t;^ (1 — 0.6m) ^ , that is, at X there is 61% Cu^, 12% Cuj, and 
27% 02x character. 



7 1-band Hamiltonians. 



The 1-band Hamiltonian for an orbital which is orthogonal to itself when displaced by a lattice 
translation, i.e. for a Wannier function, has the k-dependence of a Fourier series: 

H, (k||) = (e) — 2t (cos akx + cos aky ) + At' cos ak^ cos aky — 2t" (cos 2akx + cos 2aky ) 
+4^'-^-' (cos akx cos 2aky + cos aky cos 2akx) + 4^'-'*-' cos 2akx cos 2aky 
—2t^''^ (cos iakx + cos iaky) + ... — 2^'-®-' (cos Aak^ + cos 4aA;j,) + ... (16) 

with b eing the hopping integral between {n + l)st nearest neighbors and (e) being the average 
energy of the band. 

The conduction band shown in the right-hand panel of Fig. 3, which is given by the 8-band 
Hamiltonian for a plane dimpled by 7°, is well separated from all other bands, and this band is 
sufficiently smooth that its Fourier series converges at a reasonable pace. The results obtained 
numerically by diagonalizing Hg and Hg (with the parameters given in Ref. [18]) and subsequent 
Fourier transformation of the odd and even conduction bands are given in the first two rows of 
the following table: 
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Hg (k) / meV 
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H+ (k) / meV 


-140 
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32 
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1 „2 
16' 


1 „3 
16' 


1 „2 
16' 


1 „3 
32 



The corresponding even and odd Wannier orbitals must have fairly long range. The zero of 
energy in the table is at the average of the even and odd bands. The splitting between these bands 
will be discussed in Sect. 8. 

With decreasing dimpling, the convergence of the Fourier series (16) of the entire conduction 
band deteriorates even further. Reasons are the near-crossing of the a and 7r-bands, and the 
near-cusp of the cr-band at F. The latter can be seen in the middle panel of Fig. 3 and is caused 
by the near degeneracy of and fp. For the pure cr-band, a more useful low-energy Hamiltonian 
is therefore obtained by Fourier transformation, not of the entire anti-bonding band, but merely 
of the part near ep- We shall now give a simple, analytical derivation of such a low-energy cr- 
band Hamiltonian. The result is given in the last row of the table with r = 1/ [2 (1 -|- s)]. This 
Hamiltonian, e (k) , reproduces the a conduction band to linear order in e — ep and Eq. (22) below 
gives the second-order correction. 

Let us start with the remark, that Eq. (11) with s(e) = s{eF) is in fact a highly accurate 
1-band Hamiltonian for the a conduction band; it even describes the cusp behavior near F where 
[u,v\ ~ [1,0]. However, (11) is not an orthonormal 1-band Hamiltonian (16), because such a 
Hamiltonian must depend analytically on u and v. To obtain a 1-band Hamiltonian analytical in 
u and V, we first expand d (e) and s (e) around ep. If we only go to first order, we can solve Eq. 
(9) and obtain for the Ist-order estimate: 

£ (k) = ep+\l — d— u 

\ I + s — u 

^ £F + d-^ ^1- d- u-2rv'-' ^l + 2ru + {2ruf + ...j^ . (18) 

Note the difference between e and e. Here, d = d (ep) , d = d (ep) , s = s (ep) , r = 1/ [2 (1 -|- s)] 
and s = s {ep) . As an alternative to (12) we may use: 




In (18) we have neglected s = s (ep) and expanded 1/ (1 + s — m) for small u. Substitution of the 
cosines (6) for u and v and comparing terms with (16) yields the hopping integrals listed in the 
last row of the table given above. Here the zero of energy is taken at the average band energy 



(20) 



1/2 
1 + s 



0.32, 
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that — 
t 



and 



(21) 



For a flat plane and without s-hybridization (r=0), the 1-band Ist-order Hamiltonian £ (k) only 
has nearest neighbor interactions. For realistic s-values, hopping to second and third neighbors 
must be included, but farther hoppings can be neglected. In previous many-body calculations, t' 
was often neglected, but when it was included, t' ~ 0.3t was in fact a standard value, t" has sofar 
been neglected. 

Had we kept s/d to lowest order, we should have added to the coefficient of 2rv^ in (18) 
the term 2r(l — rf— u)s/d. Moreover, r in the table and in (20) should be substituted by r' = 

r [l - 2r(l - rf)i/rf] . 

As an example, let us take ep at the saddle-point, [u,v] = [0,1], then ep 1.53 eV, as 
found from (11), r f« 0.32 from (21), = 1 - 2r f« 0.363 from (9), rf-^ ^ 2.64 eV, from (19), 
2r (1 - d)s/d ^ 0.13 and r' = 0.28. Hence, from the table, t ^ 0.69 eV, t' /t ^ 0.28, and t" /t' ^ 
0.49. Using the first-order expression (18) to calculate the energy at ^FS, which we know from 
(11) should be 1.25 eV above the saddle-point, we find: s (^FS) -e^ = {l-d)/d= 2r'd-^ k. 1.46 
eV. This discrepancy is due to the lack of higher-order terms in (11). 

The second-order corrections to the Ist-order Hamiltonian e (k) may be obtained by expanding 
d{e) to second order and using (9) and (18) to obtain: 



e(k) 



1 [^i^)-^Ff 

2 CF - {cp + q) /2 



£(k) 



This is the exact expression for e (k) (still assuming s constant) because d (e) is a quadratic 
function and ef — (fp + frf) /2 = d/d. For the purpose of deriving an analytical 1-band 2nd-order 
Hamiltonian, we do not want to solve this equation because the result is the exact Eq. (11), 

written as: 

e(k) 




e (k) - ep 



£F - (fp + fd) /2 y £F - (fp ■ 

To obtain an analytical Hamiltonian, we must expand as follows: 



q)/2' 



e{k)-eF = [s{k)-eF] 



1-i 



£ (k) - £F 

2 CF - (fp + q) /2 



(22) 



Since €F — {ep + fd) /2 ~ 2 eV, this expansion hardly converges for e (^FS) —cf = 1.46 eV, but only 
for energies much closer to £f- This is the cusp-problem, showing up again. Our 1-band Ist-order 
Hamiltonian including 3rd-nearest-neighbor interactions is therefore appropriate for energies one 
order of magnitude closer to £f than (fp -|- e^) /2, i.e. below 200 meV. To develop an expression 
more explicit than (22) for a 1-band 2nd-order Hamiltonian which treats higher energies better 
than e (k), seems pointless considering the complication of squaring the Fourier series for e (k) — e^- 
We should also comment on the most obvious way to derive the a conduction-band Hamil- 
tonian, namely linearization and orthonormalization (3) of H (e) in (7). Using (9) and (18) it 
becomes obvious that 



n(k)-eF 



s (k) - CF 



H(eF)-eF ^ 

1-H(£f) 1 + (e(k) - £f) / (£f - £p) 



[s(k)-eF] 



1 



s (k) - CF 



Comparison with the correct expansion (22) shows that 7i (k) is merely correct to first order, like 
£ (k) , but it has longer range. In conclusion, e (k) seems to be the Ist-order 1-band Hamiltonian 
with the simplest k-dependence. 

Finally we remark that for the 8-band model, too, a 1-band Ist-order Hamiltonian may be 
derived from Ref. [18] where it is shown that the constant energy contour for an arbitrary energy 
is a polynomial of second order in x = ^ (1 — cos akx) and in j/ = ^ (1 — cos aky ) . By expansion of 
the coefficients [A(e) through 1(e)] to linear order in e — and solving, we find e — to be given 
by the ratio of two polynomials. Expansion of the denominator gives the orthonormal 1-band 
Hamiltonian. 



8 Inter-plane hopping, tj_ [k 

From the parameters of Hg and Hg given in [18] for YBa2Cu307 we derive the values of the 
integrals 

^ii = 2 i^" ~ ^^i ) ^^'^ ~ -'"2 ~ 

for hopping from orbital i in the lower plane to orbital j in the upper plane of the bilayer. We 
chose the signs of the hopping integrals in such a way that the integrals are positive if the two 
orbitals have pure cubic-harmonic character and are strongly localized. In the expression for tfj, 
the upper sign is therefore chosen for <t-<t and tt-tt inter-layer hopping integrals, and the lower 
sign is chosen for ct-tt hopping. The (tetragonally averaged and rounded) values of the inter-plane 
hopping integrals are: 



/-L /-L /-L /-L /-L /-L ^/-L /I ^/l 

SS sp pp dd zz zx,zx zy ,zy z ,zx z ,zy 

0.75 0.27 0.12 0.05 0.35 -0.30 -0.15 eV 



The negative signs presumably indicate that the corresponding hoppings mostly proceed via yt- 
trium and barium. 

For the conduction band, the most important inter-plane hoppings are the vertical tj-^ from 
Cus to Cus and tjp from Cuj in one plane to the nearest 02x or 03j, in the other plane. This is 
so, even though the Cuj character in the conduction band is relatively small. This statement may 
be verified from the left-hand side of Fig. 7, where we show 

tx(k) = i[e-(k)-e+(k)], (23) 

with the odd and even conduction bands calculated by numerical diagonalization of respectively 
and in (1), including all parameters as given in Ref. [18]. The even-odd splitting is 
seen almost to vanish along FM [M=S=(-^, ■^)] and to reach its maximum of 0.6 eV at X. This 
qualitatively follows the t;^/(l — 2rM)^ dependence of the Cuj character (14). Below, we shall 
work this out in detail. Also shown in the figure (heavy lines) are the constant-energy contour 
(FS) of the odd band passing through the saddle-point and the two neighboring ones corresponding 
to energies ±100 meV from the saddle-point. 

The right-hand side of Fig. 7 shows the same even-odd splitting and the odd-band constant- 
energy contours at, and ±100 meV from, the saddle-point, but now calculated from Hg and 
Hg . The 8-band Hamiltonian gives a more appropriate description of the even-odd splitting in 
the dimpled Cu02 bi-layer (See Figs. 1 and 2 for YBa2Cu307). The effects of the relatively 
strong hopping between 7r-orbitals on different planes are noticeable, but nevertheless, the even- 
odd splitting still attains its minimum along FM and its maximum of 0.6 eV at the saddle-point. 
That the maximum of t±_ (k) moves along with the saddle-point as it bifurcates due to dimpling 
of the plane, is caused by vertical 0^-0^ hopping. 

We now give an analytical description of the essentials of t j_ (k) for a fiat plane. Hence, we 
consider only the cr-band and neglect all inter-plane hoppings other than tj-^ and tjp. In this case. 



the Hamiltonians are (e) in (7), which depend on the mid-plane parity (_P) only through (e) . 
We therefore immediately realize that the inter-plane hopping is proporiional to v^. From (18): 

,^(k) = ( ^ )^t^iX) (^—X (24) 

eF-{ep+ed)/2 \l + s+ -u 1 + s- - u J ' \1 - 2ru J ^ ' 

where the maximum is at X [m=0,i'=1] and takes the value 

{X) = Ar' (s- - s+) — ^ 8rh ( + 2^ ) — (25) 

As usual, s = s {ep) , r = [2 {1 + s)]~ ~ t'/t, and s has been neglected. Taking again at the X- 

point of the un-split band, we get: t_L {X) k. 8-0. 322-0. 57 [jj^yIs + ^W) i 53+0 45 ^^ ~ 0.23 eV. 

Inclusion of s reduces this value by 13% to 0.20 eV. Had we instead taken ep = ^ (5^, 5^) = 2.78 
eV, we would have had s = 0, s = 0.65, and r = 0.30, so that tj_ {X) 0.17 eV. Both estimates are 
still the main part of the 0.29 eV obtained when including also those inter-plane hoppings which 
do not involve Cuj. Note that the inter-plane Cuj-Op hopping contributes as much to tj_ (k) as 
Cus-Cuj. This is because Op is a major part of the conduction-band wave function. (One might 
then wonder about the importance oitpp and t^^, but their effects tend to cancel). The fact that 
tj_ (k) in the left-hand side of Fig. 7 is not symmetric around the line XX is mainly due to the 
factor (1 - 2ru)~'^ in Eq. (24). 

Inspection of LDA calculations for other multi-layer systems reveals that the t j_ (k)-behavior 
described above is generic, as we would expect on the grounds that the Cuj (and Cu3^2_]^) orbital 
is always present and the inter-layer distance is fairly constant. What might change, are of course 
fs and tgp because they depends on the balance between the Cuj and Cu3^2_]^ characters which 
could be offset if the distance to apical oxygen became unusually short, that is, more than 0.1 A 
shorter than the 2.30 A found in YBa2Cu307. The bandstructures calculated by Novikov et al. [22] 
for bi-layer IIg-1212 and tri-layer IIg-1223, which have flat Cu02 planes separated by Ca, reveal 
that t]_ has ii^-dispersion (24) and that t]_ {X) ~ 0.25 eV, like for undimpled YBa2Cu307. Exactly 
the same is found from Mattheiss and Hamann's calculation [23] for flat bi-layer Bi2Sr2CaCu208 
(provided that one looks at the bands where they are purely Cu02-like). Finally, all calculations 
for the idealized infinite-layer material CaCu02 with fiat planes and without apical oxygen exhibit 
fc^-dispersion like 2tj_ (ky) coscA;^, exactly as expected from (24). The maximum bandwidth is 1 
eV so that, here again, tj_ {X) ~ 0.25 eV. 

A related, but less generic problem is the inter-plane hopping through insulating spacer layers, 
such as (BaO)IIg(BaO) in the Hg compounds. Here, inspection of both existing calculations 
[22, 24] for the single-layer material reveals that the fc^-dispersion is like 2tj_ (ky) coscA;^, once 
more, and that tj_ {X) 30 meV in both calculations, that is, an order of magnitude smaller 
than for hopping between multi-layers. Considering the fact, that this "backwards hopping" 
proceeds mainly through the apical-oxygen z-orbital, ii^-behavior is to be expected, because the 
04^ -orbital has m^=0 and therefore only couples to Cuj. 

t j_ (k) is the important quantity in the inter-layer pair tunnelling mechanism of P.W. Anderson 
et al. [19] for boosting T^, and ii^-dependence was in fact assumed by Chakravarty inspired by 
LDA bands. According to the inter-layer pair tunnelling theory, correlation effects block the 
hopping of single electrons between the planes, but in the superconducting state a Cooper pair 
may tunnel due to the Josephson coupling: 

Hj^-Y,Tj (k) (c^V^^-kV-ki^T + '^•c-) where Tj (k) ^ t^ (k)' /t. 

k 

BCS mean-field theory then yields the following crude estimate of the anisotropy of the supercon- 
ducting gap and T^. 

A(kf) ~ Ao(kf)+Tj(kf)/2 and ~ max Tj (k^) /4 

With our parameter values, this means that the gap for a bi-layer should be something like 
Ao (kf ) -F50meV-t;-*/ (1 - 2ruf , traced along the FS (see Fig. 7), and that ~ 300 K. 



9 Band dispersion in an anti-ferromagnetic insulator. 



It has been a challenge to understand how the valence band, or rather, its spectral density, evolves 
as a function of doping as the material changes from an anti-ferromagnetic insulator (AFI) to 
a metal, becoming superconducting at low temperature. ARPES [20] for the un-dopable AFI 
Sr2Cu02Cl2 with a single flat Cu02 plane, gave a valence-band which dispersed little from F 
(0,0) to X (tt, 0) and which rose to a nearly isotropic maximum at ^FS (f-, f-) (see Fig. 6). 
Except for the band-width, this did not agree with many-body calculations for the t — J [9] and 
t — U [25] Hubbard models, which predicted the rise to be from F to X, and the maximum to be 
flat along the AF zone-boundary (tt, 0)-(-|, -I). 

In the right-hand side of Fig. 6 we show the calculated band structure of an AFI in the simplest 
possible approximation where an external staggered field of size ±Yt^ was applied to the Cuj;2_j,2 
orbital. The results for Um = 2, A, and 6 eV are shown. In order to reproduce the experimental 
gap of 1.8 eV, one would choose m ~ 5 eV. The single-particle Hamiltonian was taken to be the 
4-band Hamiltonian [<t\ H \<t) . One observes that the calculated dispersion has exactly the same 
characteristics as the experimental one, except that the band-width is 4 times too large, as might 
have been expected. We made the points [26] that the discrepancy between the experiment and 
the many-body calculations is caused by the use of an oversimplified single-particle Hamiltonian, 
namely the 1-band nearest-neighbor Hamiltonian (the t-model), and that inclusion oft' (and t") 
should fix the problem. This has turned out to be true [27]. Let us now explain the mean-field 
results in some detail. 

For Um = (not shown), the bands are simply the anti-bonding cr-band folded into the AF 
Brillouin zone. This folding is around the t;-axis, so that u and —u become equivalent, and the 
two bands are degenerate along the t;-axis. The dispersion along the t;-axis is the 1.25 eV that 
we calculated in connection with Eq. (11) and which, on the rf-scale, is from 1 to 1 — 2r. For the 
t-model this dispersion does not exist. The relevant feature seen in Fig. 6 is that this dispersion 
along the AF zone boundary hardly changes when the staggered field is turned on. All that 
happens is that a gap of about half the size of the field (~ 50% Cu^ character) opens up while 
the dispersion essentially remains. 

The mathematics of this is really simple: The staggered field operates on Cu^ and couples 
e (k) to e (k -|- q) with q = (■^, ■^) , that is, e (m, v) is coupled to e {—u, v) . The energy dependent 
1-band cr-Hamiltonian (7) therefore becomes: 



Hie) |k) |k + q) 

(k| Q+^(l-«-T3^) fU 
(k + q| fU Q+i^fl + « 



l-|-ti + s(e) 



(26) 



Along the AF zone boundary, m = 0, we obtain: 



(2tpd)^ / v'^ \ m 



e = e, + 1 - ± -U 

e-ep \ l + s(e)J 2 

This demonstrates that the dispersion along the t;-axis is independent of the staggered field, and 
that the splitting is independent of v. Hence, the —2rv^ dispersion remains whatever it was for 
U =0. These statements are of course only true if we neglect the energy-dependences of e — fp and 
s (e) ; in reality the valence band tends to get compressed at energies closer to ep and to gain Op 
character. Solving for e yields (11) with m = and substituted by ± yt/. 

For an energy range which is narrow in comparison with the distances from ep ~ and eg, we 
may like in (18) expand the diagonal of the secular matrix [H (e) — e](e — ep) / {2tpd)^ to linear 
order around an energy at the center of interest. For consistency, we shall denote this energy ep, 
although it will usually not be the Fermi energy. In order that the resulting 2x2 Hamiltonian 
yields eigenvalues correct to linear order, the off-diagonal matrix element needs only be correct 
at ep- For simplicity, we shall neglect s, which is a reasonable approximation for the un-split and 



for the upper Hubbard band because its center is at ~4.1 eV when Um ~ b eV. The resulting 
Ist-order Hamiltonian is simply: 




(27) 



where £ (k) was given by (18) and \cd\ = (ep 



fp) (2tpd) 

|2 



is the CU(j 
0. Had we 



26f — 6p — 6d 

character (13) at ep. As previously discussed, \cdf does not depend on k when s 
gone to higher order in the e — ep expansion, the Hamiltonian would still have been (27), but with 
£ (k) and \cd\ substituted by respectively e (k) and |C(j(k)| . This form holds generally, e.g. also 
if it were derived from the 8-band model and on-site Coulomb interactions were included on other 
orbitals, when we define: 



U |c(k)p = (k)| [/ 1^ (k + q)) = J2 < (k) c, (k + q) 



From (27), the upper and lower bands are then: 



i?(k) 



e(k) + e(k + q) 



± 



'e(k)-e(k + q) 



(28) 



(29) 



If e (k) and |c (k)| are low-energy expansions, the result of (29) only applies close to ep. Hence, we 
must use different e^^'s for the upper and lower bands, as is also natural, because their characters 
are very different when m ~ 5 eV: The valence band is Op-like and the upper Hubbard band 
is CU(j-like. Remember also, that H (e) in (26) for an anti-ferromagnet describes 8 rather than 2 
bands, and that (29) can describe any of them through choice of ep and the appropriate sign. 



10 Inter-plane exchange across a bi- layer, J^. 

For anti-ferromagnetic bi-layer materials, such as YBa2Cu306, the spin-order across the bi-layer is 
anti-ferromagnetic and the inter-layer exchange-coupling constant Jj_ is an important parameter 
in many theories of high-temperature superconductivity and for understanding the origin of the so- 
called spin gap [28, 29]. From neutron scattering experiments on YBa2Cu306+j; it has been found 
[30] that, upon hole doping far into the metallic regime, the anti-ferromagnetic spin-correlations 
between the planes are more persistent that those in the planes. The value of Jj_ is under dispute, 
but it seems clear that in YBa2Cu306 the absolute value must exceed 7 meV because the optical 
spin-wave branch, which should be at energy 2^J7iyr[, could not be detected with neutron scat- 
tering up to 60 meV [31], and the intra-plane exchange constant is experimentally known to be 
-120 meV. From mid-infrared measurements of the spin-wave spectrum in YBa2Cu306 [32] it was 
recently concluded that Jj_ ~ —65 meV. SEDOR experiments on the YBa2Cu307-YBa2Cu408 
metallic compound [33], on the other hand, indicates that J_l/J|| increases strongly with decreasing 
temperature and reaches a maximum value of 0.3 just above T^. This presumably puts an upper 
bound of 50 meV on Jj_ . We shall now estimate the value of Jj_ through several implementations 
of simple mean-field theory. In all cases we assume the value Um = b eV for the self-consistent 
field applied to the Cu d orbitals. 

We consider the insulating phase, apply the staggered field ±2.5 eV inside each plane, and 
calculate the tiny difference, £ (F) — £ (AF) , in the sum of the band-structure energies according 
to whether the orientation of the staggered field is ferro- (F) or anti-ferromagnetic (AF) between 
the to planes of the bi-layer. This energy difference per Cu spin is J_l/4. Before we present an 
analytical calculation based on the results of the two previous sections, we list results of various 
numerical calculations: 



Hamiltonian: 
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TT j± j± j± j± 
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F,AF 
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0U7.2_,,2 


0U7.2_,,2 


0117.2 2 


0117.2 2 0117.2 2 

X y X y 


Jj_/meV: 


-25 


-13 


-13 


-17 


-8 -6 



These result are seen to scatter between -6 and -25 meV. Our largest calculation (LDA) used 
the charge-selfconsistent LDA potential for non-magnetic YBa2Cu306 and applied the staggered 
field to the Cuj;2_j,2 orbital in a standard multi-orbital band-structure calculation using the or- 
thonormal LMTO representation, diagonalization of the Hamiltonian, and summation over the 
occupied bands in the AF Brillouin zone to obtain £ (F) and £ (AF). Surprisingly, the same result 
( Jj_ = —13 meV) was obtained with the same brute-force procedure (F,AF), but using Hg and 
Hg for the non-magnetic part of the Hamiltonian. In the latter calculation we took advantage 
of the fact that the magnetic perturbation for inter-plane F/AF-order is even/odd so that it can 
not/only mix and H~ . In this way, the matrices to be diagonalized were only 16x16 for Hg. 
If we apply the same procedure, but start from the simplified tj-^, t^-version of H4 which includes 
inter-plane hopping via Cu s orbitals only, we find a Jj_ which is only -6 meV. The three remaining 
numerical calculations listed in the table employed a perturbative procedure labelled A(t\/U') , 
which we shall now explain. 

We use the 1-band expression (29) and calculate the change in the total energy as the neg- 
ative of the change of the total energy of merely the upper Hubbard band. This neglects 
the change of the energy of the highest, Cuj-like band. Expansion of (29) to lowest order in 
[e (k) — e (k -|- q)] /[/m |c (k) p (the result (30) will hold beyond this order) yields for the upper 
band: 

^(k).iM±£Oi±ii) + -t/|c(k)r + ^__(i^hLi^±A " 

For a bi-layer, inter-plane F-order couples (k) to (k -|- q) , while inter-plane AF-order couples 
(k) to (k -|- q). The even and odd upper bands for F-order are therefore: 

i.^(k).-"('-) + -"('- + -^) + -^|c(k)r+ ' Je^(k)-e^(k + q) - 

while the two upper bands for AF-order are 

^ ' 2 ^ 2 ' ^ ^' ^ mU\c(k)f V 2 

The difference in total energy per Cu spin thus works out as: 

-i J, ^SiF)-S (AF) = i (E^ (k) + E'^ (k) - E+ (k) - i?" (k)) = ('-±^1^^^) , 
^ z \ fjiU \c (k) I / 

where () denotes the average over the Brillouin zone. As usual (23), t±_ (k) is half the splitting 
between the odd and even bands for U=0. Expression (30) in real space is, in fact, the familiar 
one for Anderson super-exchange. 

The results labelled 4:{tl/U) were now obtained by diagonalizing the appropriate non-magnetic 
4x4 or 8x8 Hamiltonians (k) and H~ (k) numerically, extracting the odd-even splitting to 
form t±_ (k) , the eigenvectors to form rnU |c(k)| , and summing over the Brillouin-zone (30). The 
result obtained in this way for the simplified 4-band Hamiltonian ( Jj_ =-8 meV) is in accord with, 
but slightly larger than that obtained with the direct procedure for the same Hamiltonian. The 
result obtained for the full 4-band Hamiltonian (Jj_ =-17 meV) is twice as large. Finally, since 
the conduction band of the 8-band model has Cu^;^ and Cuj,^ components besides Cuj;2_j,2, we 
did a calculation in which |c(k)p denoted the sum of the characters for all three Cu rf-orbitals. 



rather than just that of Cuj;2_j,2. The corresponding decrease of the denominator in (30) raised 
J_L to -25 meV. 

We thus conclude that, for the purpose of calculating Jj_ in the mean-field approximation, Hg 
simulates the full LDA Hamiltonian very well, the perturbation-expression (30) for Jj_ is accurate, 
the assumption of inter-plane hopping via only Cuj under-estimates Jj_ by a factor 2, by-mixing 
of Cu diT character enhances Jj_, and our best estimate of Jj_ is ~-20 meV. 

With expression (30) we may evaluate Jj_ analytically for H4, using (13) for the Cuj;2_j,2 
character and the Ist-order expressions (24) and (25) for t± (k) in the t^^,, ^^-approximation. We 
obtain: 

-J. = (r)' f 1 + 2 ^ + . 4|^I^ . 0.14 = 6 meV (31) 

mU\cdf V 9 J 5 eV-0.57 ^ ' 

where t±_ (X) was given in (25) and the last factors come from the Brillouin-zone average: 

{(TriW)) = (''M'-'-''-J) = (l)^(-r-> 

The numerical estimate for — Jj_ was obtained by using the expansion energy ep = e (^, ^) = 2.78 
eV, which is at the center of the U =0-band, and for which we can safely neglect s. The -6 meV 
is in good agreement with the values obtained numerically for the same Hamiltonian. Since the 
ii'^-dependence dominates the ru-dependence in the average over the Brillouin-zone, expression 
(31) for Jj_ is fairly insensitive to details of the dispersion. If for t±_ (X) we merely insert the LDA 
value of 0.25 eV, expression (31) in fact yields: Jj_ —13 meV, which is the LDA result. 

Barriquand and Sawatzky [34] recently obtained Jj_ ~ —56 meV from a real-space evaluation 
of (30). One reason for the discrepancy with our value is, that these authors only considered the 
perpendicular hopping integral t^^ and took it to be 2.5 times larger than our value (as given in 
Sect. 8). It seems to us, that had the perpendicular hopping been dominated by tpp, then t±_ (k) 
would have been nearly independent of k, in contradiction with the even-odd splittings found in 
all LDA calculations. 

In a numerical (F,AF)-calculation, we can simulate hole-doping crudely by moving the Fermi 
level into the valence band. At the same time, the mean field should be reduced in such a way that 
it vanishes at a hole doping of about 0.2 per plane. Fig. 8 shows the result of such a calculation 
of Jj_ for many different values of the field. The persistence of the inter-plane anti-ferromagnetic 
coupling with increasing doping is clearly seen, until a point where Jj_ changes sign and the 
coupling becomes ferro-magnetic. 
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FIGURE CAPTIONS 

Fig. 1. LDA energy bands, orbital projected bands, and Fermi surface for YBa2Cu307 and 
k,=0 [3]. 

Fig. 2. Same as Fig. 1, but for k;, = j. 

Fig. 3. Specification of the 8-band Hamiltonian and synthesis of its band structure. All 
energies are in eV. See text. 

Fig. 4. 02j;-Cuj;2_j,2-03j, conduction-band states at k = (-^,0) and (5^,5^) , schematic. 

Fig. 5. LDA wave-functions ^surfaces of constant |i/'(k,r)p^ for the odd plane conduction 

band in YBasCuaOr. Left: k = (f , 0, O) X. Right: k = (^, ^, O) iPS [21]. 



Fig. 6. Left: Non-magnetic and anti-ferromagnetic Brillouin zones and relation between the 
(kx,ky) and [u,v] coordinate systems. Right: Energy bands from the 4-band Hamiltonian in the 
presence of anti-ferromagnetically staggered fields acting on the Cuj;2_j,2 orbital and of strengths 
±1, ±2, and ±3 eV. 

Fig. 7. Contours of constant t±_ (k) = ^ [e~ (k) — (k)] calculated for the 4- and 8-band 
models. The heavy lines give the odd constant energy surfaces passing through the saddle-points, 
as well as the two neighboring contours corresponding to the energies ±100 meV from the saddle- 
point. 

Fig. 8. Jj_ as a function of hole doping calculated for many different values of the staggered 
field. The single-particle Hamiltonian was H4 tj-g,tjp and the procedure was (F,AF) as explained 
in the text. 
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